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Abstract 

We study a nonlinear system of partial differential equations in which a complex field (the 
Higgs field) evolves according to a nonlinear Schrodinger equation, coupled to an electromagnetic 
field whose time evolution is determined by a Chern-Simons term in the action. In two space 
dimensions, the Chern-Simons dynamics is a Galileo invariant evolution for A, which is an 
interesting alternative to the Lorentz invariant Maxwell evolution, and is finding increasing 
numbers of applications in two dimensional condensed matter field theory. The system we 
study, introduced by Manton, is a special case (for constant external magnetic field, and a point 
interaction) of the effective field theory of Zhang, Hansson and Kivelson arising in studies of the 
fractional quantum Hall effect. From the mathematical perspective the system is a natural gauge 
invariant generalization of the nonlinear Schrodinger equation, which is also Galileo invariant 
and admits a self-dual structure with a resulting large space of topological solitons (the moduli 
space of self-dual Ginzburg-Landau vortices). We prove a theorem describing the adiabatic 
approximation of this system by a Hamiltonian system on the moduli space. The approximation 
holds for values of the Higgs self-coupling constant A close to the self-dual (Bogomolny) value of 
1. The viability of the approximation scheme depends upon the fact that self-dual vortices form 
a symplectic submanifold of the phase space (modulo gauge invariance). The theorem provides 
a rigorous description of slow vortex dynamics in the near self-dual limit. 

1 Introduction and statement of results 

In this article we study vortex dynamics in a nonlinear system of evolution equations (jl.5p intro- 
duced by Manton (1997). This system is in fact a special case of an effective field theory for the 
fractional quantum Hall effect (the Zhang-Hansson-Kivelson, or ZHK, model). In addition it is 
a natural gauge invariant generalization of the nonlinear Schrodinger equation, possessing impor- 
tant structural features (Galileo invariance and self-dual structure with existence of related moduli 
spaces of solitons) which make it interesting to study for mathematical reasons. After introducing 
the system under study, and putting it into mathematical and physical context, we explain the 
necessary background material in order to state our results, which appear in Q1.71 

1.1 Chern-Simons vortex dynamics 

We start by motivating the study of Manton's system from the mathematical perspective, before 
going on to show that it is equivalent to a special case of the ZHK model, and discussing its physical 
significance. 



1.1.1 Manton's system on R 2 : mathematical context 

To introduce Manton's system, we start with the nonlinear Schrodinger equation on R 2 : 

i^- = -A$- -(1- |$| 2 )$, (1.1) 
dt 2 y 1 1 ; ' K J 

to be solved for $ : R x R 2 — > C; A is a positive number. This has the following properties: 

(i) it defines a globally well-posed Cauchy problem, 

(ii) it admits topological soliton solutions, the Ginzburg-Landau vortices, and 

(iii) it is invariant under the group of Galilean transformations. 

Manton's system is a generalization of (|l.lj) . sharing these properties, which describes the 
evolution of a complex field $, coupled to a dynamically evolving electromagnetic potential A = 
AqM + A\dx l + A2dx 2 . On R 2 the system reads explicitly (writing (a, b) = ^kdb): 
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In addition to (i)-(iii) above, this system has the following mathematical properties: 

(iv) it is gauge invariant, 

(v) self-dual structure and a large space of topological solitons (see §1 .6j) . 

These properties make the study of vortex dynamics in Manton's system interesting, since the 
self-dual structure makes a rigorous analysis possible when the vortices are arbitrarily close (see 
§1.6II1.7|1 . The proof of our results makes use of special mathematical features present due to self- 
duality which are explained in £j3j these features include complex and symplectic structures on the 
soliton moduli space, and a foliation of the phase space which we call the Bogomolny foliation. 

1.1.2 Equivalence of Manton's system and a special case of the ZHK model 
The system (II. 2j) can be derived from the action S = cq f s(A, &)d 2 xdt, where 

s(A, $) = -e^A^Ap + (d t - iA )$) +A + W k d 3 A k \ 2 + |(S,- - iA^\ 2 + ^(1 - |$| 2 ) 2 

where Greek indices run over {0, 1,2} for space-time tensorial quantities, Roman indices run over 
{1,2}, and e ,J ' up , e J are the completely anti-symmetric symbols and the summation convention is 
understood. This action is one of a class involving the Chern-Simons term e^ up A^d u A p , see [17\ 122] 
for a review. It is characteristic of these theories that variation of the action with respect to Aq 
gives a constraint equation involving the magnetic field, in this case the final equation of (|1.2|1 . 
This equation is analogous to the Gauss law in ordinary Maxwell theory, and is referred to as 
a constraint because the previous (dynamical) equations in (|1.2|) imply that its time derivative 
vanishes (exactly as do the dynamical Maxwell equations for the Gauss law). This constraint 
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means that many apparently different actions give rise to the same Euler-Lagrange equations: in 
particular we can replace the above action density with 



s(A, $) = -e^A^Ap + (**, (d t - L4 )$) + A + |(V - zA)c|>| 2 + — j-(l - |<D| 2 ) 2 . 

We now introduce the ZHK action S zllK (a, <£; A ext ) = cf s ZUK d 2 xdt and show that Manton's 
action S is in fact a special case of S zllK ; essentially the same observation appears also in [22, page 
54] and [381 Section 4.8]. The ZHK action is the action for a mean field description of the quantum 
Hall effect. This effect refers to the current Jj = aj/,E^ xt produced in an effectively two dimensional 
system of electrons in a strong transverse magnetic field, by application of an applied electric field 
E^ xt . In the right experimental situation the conductivity tensor is found to be off-diagonal 
(i.e. o\\ = = (J22), with the non-zero entries a 12 = — 021 = fe 2 /h, where / is an integer, or a 
fraction, for (respectively) the integer and fractional quantum Hall effect. This quantization of the 
values of a 12 means that as the number of charge carriers is increased there is no corresponding 
increase in the current - it lies on a plateau - at least until the number of carriers is sufficiently 
greatly increased, at which point the conductivity moves to another of the quantized values, and 
the current moves to another plateau. In the mean field description the field $ interacts with an 
external (applied) electromagnetic potential A ext and a "statistical" potential a, according to: 

Szhk = ^ afl d u a p + (d t - ia ~ i^)*> + ^-|(V - ta- iA ext )<S>\ 2 
2 2m 

+ J (I - \<S>(x)\ 2 )V(x - x')(l - \$(x')\ 2 )d 2 x', 

(see [52], [T71 Section 4.6], or [53], Equations (7)-(8)], taking note of the published erratum for the 
latter reference). To reduce this to s we consider the case of a constant external magnetic field 
B ext = diA^ 1 — dzAf* with A^f 1 = 0. (The standard configuration in quantum Hall experiments 
involves a strong transverse magnetic field applied to an effectively two dimensional electron gas, 
with relatively small electric potentials applied along one of the planar directions.) Define A = 
a + A ext . Now check that 

e flvp a fJ ,d v a p = a (<9ia 2 - d 2 ai) - ai(d t a 2 - d 2 a ) + a 2 (d t ai - c^ao) 
= A)(diA 2 - d 2 A 1 - B ext ) - (A\ - Af t )(d t A 2 - d 2 A ) 

+ (A2-Al xt )(d t A 1 -d 1 A ) 
= e^A^Ap - 2A B ext + d t (Af t A 2 - A e 2 xt A{) + 3 x (A Q Af *) - d 2 (A Af *) 
= e^A^Ap - 2A B ext + ^d^Af* A p ) 

and deduce that s ZHK — s ZHK is a derivative, where 

Szhk = ^PA^Ap - KB ext A + (dt - iAo)*) + ^|(V - iA)$\ 2 

+ J (l-\$(x)\ 2 )V(x-x')(l-\$(x')\ 2 )d 2 x'. (1.3) 

Now recall that derivatives in the action density do not affect the corresponding Euler-Lagrange 
equations (they are null Lagrangians) . It follows by comparing s and s ZHK that the equations of 
motion for S^hk will be identical to those of Manton if we choose V(x) = (A + l)<5(x)/4, k = 
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—2,B ext = 1/2 and m = 1/2. Therefore, we conclude that, at least as far as the classical equations 
of motion are concerned, the ZHK model with these values is the same as Manton's system in the 
case of 

• a constant external magnetic field of appropriate value, and 

• a point interaction V(x) oc 5(x). 

We now discuss the physical interpretation of the model in the fractional quantum Hall context. 
There is a microscopic model, due to Laughlin, which explains the observed phenomena in a well- 
accepted way in terms of a new phase of the two dimensional electron gas (for low temperature 
and high magnetic fields), with ground state described by the Laughlin wave function. There is 
an energy gap in the spectrum, so that the excitations above this new ground state have strictly 
positive energy - the Laughlin quasi-particles and quasi-holes, which have fractional statistics and 
fractional charge. It is this fractional charge, combined with the explanation of the integer quantum 
Hall effect, which gives rise to the fractional quantum Hall effect. The effective field theory proposed 
in [53], and reviewed at length in [52], gives a mean field description which is not expected to be 
accurate on microscopic length scales, but which does give an alternative explanation of all the 
main observed phenomena. In this mean field theory, the elementary excitations are described by 
the topological vortices, which are endowed with the same fractional charge and fractional statistics 
as the Laughlin quasi-particles. (The Chern-Simons term for the statistical gauge field a in the 
action serves to change the statistics in the well known way explained in [521 [2] ) . It is understood 
that, in the mean field picture, it is the pinning of vortices which explains the observed plateaus 
in the Hall conductance ( [MJ [43J [52j ) . Thus a good understanding of vortex dynamics in the ZHK 
model should be useful to gain a better explanation of the phenomena within the context of the 
mean field approach. Needless to say there is still much work to be done to go from the results 
of this paper, to results which would apply directly to the experimental situation: even apart 
from issues like the spatial domain and the real values of the coefficients in the model, it will be 
necessary to treat the applied electric potential which produces the Hall current flow. This means 
that in the above derivation we should allow for an external electric field E? = —d^Ao, in addition 
to the static magnetic field, and investigate its effect on the vortex motion. Since the magnetic 
field is very strong in the experimental situation, it should be reasonable to treat the electric field 
perturbatively. 

To conclude the discussion of the motivation for our work, the system (jl.2p is one of a class of 
dynamical Chern-Simons vortex models whose study is mathematically interesting (due to proper- 
ties (i)-(v) above), and which is physically relevant (as we have just discussed). The use of such 
models in condensed matter applications is phenomenological, so the precise Lagrangian and many 
values of the coupling constants, etc. are not precisely known. (Actually, in [53^ 152] the ZHK action 
(jl.3p is derived formally from an ostensibly microscopic, second quantized, action. However, this 
microscopic action itself seems to have a phenomenological character, since it involves excitations 
which are not fundamental electrons, but rather collective excitations - see the dicussion following 
(2.6) in [43]). In any case our main result does provide a rigorous basis for understanding vortex 
dynamics in a prototype for a class of theories which are of interest in two dimensional condensed 
matter theory. The adiabatic limit system (|1.22j) which we derive for the vortex dynamics cannot 
usually be written down explicitly, but as discussed in remark 11.7.41 the behaviour of some of its 
solutions can be understood in many cases, and thus information on the dynamics of vortices can 
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be deduced within the framework of this approximation. There are reasons to hope that qualitative 
features of the motion in this limiting situation will have a wider validity: see remark fl.7.41 

As a final comment on the quantum Hall effect, there is another type of soliton - a nonlocal 
Skyrmion - which appears in treatments of the ferromagnetic properties of quantum Hall samples 
(see [121 [501 [18] and [H] for some analytic properties of these Skyrmions in a particular case) . 

1.1.3 General physical context for Chern-Simons models 

There has been a fairly long standing interest in systems of the type (jl.2p in the physics literature; 
we give a brief summary and refer the reader to |23[ \T7\ [22] for detailed reviews. The study of 
Chern-Simons dynamics in 2+1 dimensional Maxwell and (non-abelian gauge) theories was started 
in the early 80's (see e.g. |11| ) and the incorporation of vortices into this dynamics (in systems 
with coupling to a nonlinear Schrodinger equation) has been studied since at least the early 90's 
by theoretical physicists (see the review |23] for early work on Chern-Simons vortices). The reason 
for this interest is both because (i) the Chern-Simons models are used widely in condensed matter 
physics in descriptions of the quantum Hall effect and high T superconductivity, and (ii) because 
they provide a useful scenario in which to probe certain complex issues in field theories. 

Regarding the first point, there are various time-dependent models for magnetic vortices but 
at very low temperatures it is argued ([3j [44]) that the motion should be non-dissipative so the 
usual Eliashberg-Gorkov equation is not appropriate, and the Chern-Simons coupled to Schrodinger 
vortex dynamics is widely used instead in the condensed matter literature, both in superconduc- 
tivity and the quantum Hall effect; see [331 Section 10.7], [32l Chapter 6] for general discussions, 
in addition to the references for the ZHK model in the previous section. (Relativistic invariance 
is broken in these condensed matter applications, so the corresponding relativistic abelian Higgs 
model, whose vortex dynamics are studied in [45J, is not appropriate. The main application which 
has been suggested for the relativistic dynamics appears to be cosmic string evolution.) There have 
been explanations offered for the wide occurence of Cherns-Simons types models in two dimensional 
condensed matter applications in terms of universality features of large scale effective actions for 
two dimensional interacting electronic and magnetic systems with spin ( [16\ Section 3]). 

Regarding the second reason for interest in these models, it was realized in the 1980's that in 
two dimensions there were possible quantum statistics other than the usual fermionic and bosonic 
types - anyons, are two dimensional quantum particles undergoing an arbitrary phase shift on 
interchange. Furthermore, composite objects made up from charged particles orbiting vortices (or 
flux tubes) have fractional spin and statistics ([51]). In [15] the authors study the quantum theory 
of a Lagrangian which is closely related to (j!.12[) . and use it to investigate the quantization of 
solitons, quantum statistics and anyons in a rigorous quantum field theory setting. 

1.2 Organization of the article 

The article is organized as follows. Our main aim is the study of vortex dynamics in the Chern- 
Simons-Schrodinger system with spatial domain a Riemann surface, so we start in the next section 
by writing down the equations in this case, and then giving necessary background including a 
discussion of the self-dual vortices in §1.61 We then state our main result, theorem 11.7.21 which 
describes the adiabatic approximation of vortex motion in the self-dual limit. This is proved in $2] 
following a strategy explained in the context of a simple model problem in £11.81 The proof uses 
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some specialized identities related to the self-dual (or Bogomolny) structure, presented in Sj3] (which 
may be read separately). Various subsidiary facts and lemmas are given in the appendix. 



1.3 The equations on a surface 

The dependent variables are a complex field &(t, x), and an electromagnetic potential 1-form 

A dt + Aidx 1 + A 2 dx 2 . 
This 1-form determines a covariant derivative operator 

D = (D , Di,D 2 ) = (g- t - iA , D U D 2 ) = (— - iA , Vi - iA lt V 2 ~ %A 2 ), (1.4) 

which in turn determines the electric field E = Ejdx J and magnetic field B(t,x) via (jl.6p : all these 
fields are defined for (t,x) € K x S where S is a two dimensional spatial domain, taken to be 
a Riemann surface with metric gjkdx 3 dx k , area form dfi g and complex structure J : T*S — ► T*S 
(where j, fc, . . . take values in {1, 2} and we use the summation convention). Introducing a covariant 
Laplacian operator by 

-A A d> = -—^Djtfiy/d^jDiQ) 
Vdet g 

(using a local frame and coordinates), the equations are 



- L4 )$ = -A A cD _ ^(i _ |$|2)$ (1.5) 
B= 1 -(l-\$\ 2 ). 

The electric and magnetic field can be combined to give the space-time electromagnetic field 

F^dx 11 A dx v = Ejdt A dx j + Bd\i g . 

This two form is obtained as the commutator of the space-time covariant derivative (|1.4p which 
mediates the coupling in (jl.5p : 



[L> M , D v }<$> = -iF^Q, where F ofc = E k , and ^F jk dx j dx k = Bd^ g . (1.6) 

(Greek indices run through 0, 1, 2 and Latin indices through 1, 2 only. Bold face is used to indicate 
the spatial part of a vector or one- form etc., except in Jj3] where time does not appear at all.) 

We now describe this set-up briefly in geometrical terms. Assume given a one dimensional 
complex vector bundle L — > S, with a real inner product h locally of the form (a, b) = h^tab, and 
corresponding norm |a| 2 = (a, a); if we employ a unitary frame over some chart then (a, b) = ^tab. 
We are then solving for an S 1 connection on the bundle LElxi-»KxS, with associated 
covariant derivative D, and a section <3? of L. To be more explicit, fix a smooth connection on L 
determined by a covariant derivative operator V, so that the spatial part of D, which will be written 
D, takes the form Dj = Vj — iAj for a real 1-form A = Ajdx^ £ il^(S); here V is independent of 
time. (It is generally not possible to choose V to be flat, and it will have a curvature, determined 
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by a function b such that [Vj, V k]&dxi dx k = —ibdfi g ^; it is always possible to choose b = const. , 
and we will do this throughout.) In any case, with this procedure the space of connections on L 
can be identified with the space of real one-forms. Then at each time t G K we are solving for a 
section $(£) of L, a 1-form A(i) = J 4i(i)(ix 1 + ^(^cfcc 2 on E, and a real valued function Ao(t) on 
E. The electric field is given by 

E= dA l _dAo 
3 dt dxi 

and the magnetic field by 

Bd[i g = bdfig + dA. 

(Here, and elsewhere, we write d in bold face when it is necessary to indicate that only the spatial 
part is taken.) The 2-form —iEjdt A dx 3 — iBd/j, g is the curvature associated to the space-time 
covariant derivative D, as in (jl.6p . For the case E = M 2 , the system was proposed by Manton 
(1997), who derived it as the Euler-Lagrange equation for the Lagrangian (|1.12j) . 

Notation 1.3.1 We shall always consider conformal co-ordinate systems on E in which the metric 
is of the form g = e 2p ((dx 1 ) 2 + (dx 2 ) 2 ) and the volume element is then e 2p dx 1 A dx 2 . On functions 
the Hodge operator acts as *f = fdfj, g = fe 2p dx 1 A dx 2 and * 2 = 1, so that *du = e _2p (|^f — f^-) 
for 1-forms uj. On 1-forms *{u\dx 1 + u^dx 2 ) = io\dx 2 — ^dx 1 , which is just the negative of the 
complex structure J, represented in conformal co-ordinates by the anti-symmetric tensor J\ with 
j\ = —1, J 2 = +1, the other components being zero. Correspondingly we decompose a one-form as 
u> = ui^'^dz + uj^ ,l ^dz; in particular for the derivative df = dfdz + dfdz, with Bf = + 
and 

D$ = D (1 ' 0) $ + D 1 - ' 1 ^ = d A $dz + 5 A <f>dz, 

with 8a& = \ ((Vi — iA\) +i(V2 — 2^-2)) 3* etc.; see 53 For a 1-form A we write the co- differential 

d*A = —divA, with divA = e 2p (^r + an d the Laplacian on real functions is Af = 

e~ 2p Q^ig x i , (with the summation convention), and on sections of L the covariant Laplacian is 
-Aa$ = e~ 2p (D 2 + when a unitary frame is used. The operators div,*d,A (resp. all 
depend on g (resp. g,h), but this is not indicated as g,h are fixed, and similarly dependence of 
constants in estimates on (E, g) and h will be suppressed throughout the article. 

Notation 1.3.2 We are dealing with sections of smooth vector bundles V over E with an inner 
product (•,•) induced from the Riemannian metric g and the metric h on L in the standard way; 
since g, h are fixed throughout they will not be indicated. Thus, for example, 

|D$| 2 = ((DiS, .Di$) + (D 2 $,£>2$))- 

We write £l°(V) for the smooth sections of V and Q P (V) for the smooth p-forms taking values in 
V. We will make use of the Sobolev spaces H S (V) of sections ofV whose coefficient functions (in 
any frame over any open set £1 C Ej lie in the standard Sobolev space H s (£l); the corresponding 
Sobolev space of V -valued p-forms is denoted H s (£l p (V)). In JF] and we shall generally omit 
explicit reference to the vector bundle, since this is usually clear, and write H s in place of H S (V) 
etc. (and \\ ■ \\jjs for the corresponding norms). However if it is necessary to emphasize that time 
is fixed, and the norm is taken over S ; we shall write H S (T,). 

Further notational conventions are given in the appendix and in particularly in relation to 
the complex structure (see also the textbook [24, §9.1] for a treatment of the background material). 
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1.4 Existence theory for the Cauchy problem 

Inherent to the system (jl.5p is the property of gauge invariance: let x(A x ) be a smooth real valued 
function, then (A, <£) is a smooth solution if and only if (d\ + A, &e lx ) is. This introduces a large 
degeneracy to the solution space which may be removed by a choice of gauge in various ways. We 
will adopt here the following gauge condition which involves the time derivatives A, <3?, of A, <!>: 

div A - (£$, <£>} = e~ 2p (diii + <9 2 i 2 ) -(**,*) = 0. (1.7) 

We make this choice because it allows a convenient description of the complex and symplectic 
structures on the moduli space of vortices (see remark 11.6.31 and $3]), and also is useful in the 
derivation of energy estimates for the time derivatives (see §2.21 and §2.3h . In this gauge global 
existence can be stated as follows: 

Theorem 1.4.1 (Global existence in gauge (jl.7p ) Consider the Cauchy problem for fjl .5j) with 
initial data $(0) G 7f 2 (X) and A(0) G -fT 1 (S). There exists a global solution satisfying (ll.7p and 
the estimate 

mt)\H*m < ce ae0t (1.8) 

for some positive constants c,a,/3 depending only on (E,g), the equations, and the initial data. The 
solution has regularity $ G C([0, oo); H 2 (Y,)) n C 1 ([0, oo); L 2 (£)) and A G C 1 ([0, oo); ^(E)) . 7/ 
i/ie initial data are smooth, then the solution is also smooth. 

It is explained in appendix IA.3I how to derive this theorem from the global existence result of 
[13] . which is stated in another gauge. Bounds of the type (|1.8p were derived in |T0] for the cubic 
nonlinear Schrodinger equation on R 2 , by means of the inequality 

|n| L oo <C[l + v/ln(l + \\u\\ H 2 )], (1.9) 

valid for u G H 2 (R 2 ) and with C = C(||u||#i). The proof of global regularity for (|1.5p depends 
on a covariant version of this inequality (given in lemma IA.11|) , and a careful treatment of various 
commutator terms \D^,D V \ which indicates that they have a comparable strength to the cubic 
nonlinear term. 

In conclusion, theorem 11.4.11 provides a global solution which is a continuous curve in the space 
Ti.2 where for s G K we define 

H s = {(A, $) G H'-^T.) x 77 S (S)}, (1.10) 

with the corresponding norm || • \\n s . From now on we will consider only (A, <3?) which lie (at a 
given time) in the space Ti.2- The gauge group at fixed time is given by 

G^igeH^^S 1 )} (1.11) 

and acts on H2 according to g • (A, 3>) = (A + g^dg, &g). (Restricting to the set where $ is not 
identically zero the action is free and gives a principal <y— bundle structure. The gauge condition 
(ll.7p can be then regarded as giving a connection - i.e. a family of horizontal subspaces - on this 
bundle.) 
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1.5 Variational and Hamiltonian formulation 



The equations (|1.5p can be derived formally as the Euler-Lagrange equations associated to the 
functional 

S(A,$) = 1- [ -AAF+((i$,DQ$)+A + 2v x (A,$))dtdfj, g , (1.12) 

where 

v x (A, $) = i(i? 2 + |D<D| 2 + ^(1 - |<I>| 2 ) 2 ) (1.13) 

is the density of the Ginzburg-Landau static energy. (The parameter A is a positive real numbers). 
Although S is not manifestly gauge invariant it changes by an exact form under gauge transforma- 
tion, and the Euler-Lagrange equations (|1.5|) are gauge invariant. Vortices are critical points of the 
static energy 

Va(A,<&) = / v x (A,$)d» g , 

as will be discussed further in the next section. 

To see that the system (jl.5p is Hamiltonian, observe that there is a complex structure on the 
phase space Ti.2 given by JT : (A,$) = (— JA,i<&) which allows the introduction of a symplectic 
structure Q(u, w) = (3v, w) where ( • , • ) is the L 2 inner product. Using this symplectic form the 
system (jl.5p . in temporal gauge Aq = 0, is a Hamiltonian flow generated by the Hamiltonian 
functional V>(A, $), which was just defined. (A short calculation reveals that the third equation 
of (|1.5|) is preserved by the evolution, and as such is really only a condition on the initial data. It 
will be referred to as the constraint equation.) 

1.6 Self-dual vortices and dynamics in the limit A — > 1. 

The system (11.51) admits soliton solutions, called abelian Higgs, or Ginzburg-Landau, vortices, 
which are energy minimizing critical points of the static energy functional Va(A,$). We now 
discuss these solutions and their uses in understanding the dynamical system (I1.5P via the adiabatic 
approximation. There is a special case, A = 1, in which the adiabatic approximation is particularly 
powerful because the space of vortices is then unusually large - large enough that the motion on it 
can provide information on the dynamical interaction of several vortices. We call this the self-dual, 
or Bogomolny, case, and the corresponding solutions are called self-dual vortices. Now for such a 
solution, (A, <£), with a given value of the topological integer N, (the degree of L), the field $ will 
have N zeros, counted with multiplicity. Each of these zeros can be thought of as the centre of a 
vortex. Thus the static solitons can be thought of as a nonlinear superposition of N vortices which 
do not interact. This was first fully understood in the case that X is the upper half plane with 
canonical metric, when the equations were solved exactly by Witten (1977) by reducing them to the 
Liouville equation. In general it is still possible to make a reduction to a nonlinear elliptic equation 
of Kazdan- Warner type, whose solutions can be completely parametrized although not explicitly 
given. Following this, Taubes proved an existence theorem when X is the Euclidean plane (Jaffe 
and Taubes 1982), and Bradlow (1988) did likewise for £ a compact Riemann surface, proving the 
following: 

Theorem 1.6.1 (Existence of vortices on a surface, |8j) If the area of a closed Riemann sur- 
face |S| is such that |S| > 4ttN the Bogomolny bound is saturated: in fact the minimum value irN 
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of V, where 



V : H 2 -> R (1.14) 
V(A, $) = Vi(A, $) = ±Ub 2 + |D<D| 2 + ^(1 - |<D| 2 ) 2 ) c^, 

is achieved on a set Sjy C °/ pairs (A, <£) which solve the Bogomolny, or self-dual vortex, 

equations: 

d A $ = o, £-I(i-|$| 2 ) = o. 

These minimizers will be referred to as the self-dual vortices, or just vortices. The quotient of Sn 
by the gauge group Q can be identified with Sym N (T,), the symmetric N-fold product o/E, via the 
mapping which takes to the set of its zeros. 

Remark 1.6.2 (Interaction and stability of vortices) The physical interpretation of theorem 
11.6. H is that for A = 1 the vortices do not interact; see [28] for a discussion of this, and some related 
conjectures, and [19J for some stability theorems. 

Remark 1.6.3 (Bogomolny structure and Bogomolny operator) The structural feature of 
V which makes theorem 11.6.11 possible was identified by Bogomolny in [7]. In this instance it 
amounts to the fact that if we introduce the Bogomolny operator B to be the nonlinear operator 
which maps (A, $) i-> (B — |(1 — |<3?| 2 ), 8a$) then 

V = \j \B(A,^)\ 2 d^ + nN 

(see $3] for more information in this regard) . Also see [8] for higher dimensional versions of this 
decomposition, and [21] for generalizations to solutions with non-vanishing electric field. 

Remark 1.6.4 (Geometry of moduli space) Quotient spaces of the type arising in theorem 
11.6.11 are usually known as moduli spaces: in this case we define the moduli space A4n to be the 
space of gauge equivalence classes of self-dual vortices, so that A4n = Sym Ar (S). We call the 
space Sn the vortex space and proj : Sn — ► A4n the natural projection which takes (A, <£>) to 
its gauge equivalence class [(A, <&)]. The space M.n inherits both a metric (induced from the L 2 
metric) and a symplectic structure and is a Kaehler manifold (see [9]). Explicitly, we can identify 
the tangent space to Mn with solutions (A,<3?) of the linearized Bogomolny equations which also 
satisfy the condition (|1.7j) . The complex structure and symplectic structure on M.^ are then given 
by restricting the formulas given in the previous section to such (A,$), and consequently we will 
use the same notation, J and O, for these objects. The existence of this complex structure on 
M m can be seen very clearly in the formulas in £j3j in which complex notation is used to combine 
the linearized Bogomolny equations with (II. 7h into a manifestly complex linear operator T>^, for 
tp = (A, <I>) G Sn- This can all be summarized by saying that we have an identification 

T[ip]M.N ~ KerV^ = {(A, $) : DB^[A, = 0, and flTTD holds}. (1.15) 



10 



1.7 Statement of the adiabatic limit theorem 



In order to define the adiabatic limit system, we now define a Hamiltonian function M.n K 
by restricting the energy Va to the space of vortices, and observing that by gauge invariance this 
actually gives a smooth function on the quotient space Mn- The corresponding Hamiltonian flow 
determines the slow motion of vortices for A close to 1: 

For e = | A — 1| sufficiently small, the system f 1 1 . 5 j) can be approximated, for times of 
order by the Hamiltonian flow on the phase space Mn = Sym N (Yi) associated to the 
Hamiltonian function VaIa^jv the symplectic form fl 

We now move towards a precise formulation of this in theorem 11.7.21 Since we are interested in the 
regime in which |A — 1| <C 1 it is useful to introduce a large parameter 



and let also, for A ^ 1, 



A - 1 , 
" = iI3T|==" (1.17) 



< •;■>;■>< M \ ; . \ \ r j .cm <iu • I i mi • i )\ / 

leading to the following rescaled equations: 



(also defining a = for A = 1 where necessary). We rescale time by r = ^, and Aq similarly, 



^=^(- ft S + (i *,J> 1 *)) + ^, (1.18) 

i(— - iA Q )$ = jx(-A A * - ^(1 " l^l 2 )^) " \ (1 " l^l 2 )^- 

It is also natural to separate the energy Va into the (main) self-dual piece V = Vi, and a perturbation 
term proportional A — 1. Under the rescaling just introduced, the energy rescales by a factor \x, 
leading us to consider the Hamiltonian H = + U, where V = Vi is as in (|1.14p . and the energy 
correction away from the self-dual, or Bogomolny, regime is given by 

LT(*) = ! Jjl-m 2 ) 2 dv 9 . (1.19) 

The rescaled equations (|1 . 18[) can be written as a Hamiltonian evolution for ip = (A, $) in the 
form 

JJ^ = + U' + S(dA ,iA <S>) (1.20) 

OT 



where J is the complex structure introduced at the end of §1.51 

I^Aidx 1 + A 2 dx 2 , $) = {-A 2 dx l + Axdx 2 , i$) (1.21) 

with A = 3£ * = |£. 
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Remark 1.7.1 (Explicit formulation of adiabatic limit system) We now write the equations 
for the adiabatic limit system in an explicit way which will be useful later. The function U is clearly 
gauge invariant and defines by restriction a smooth function u on Mn- Now recall ()1.15f) : under 
this identification, the gradient of the function u on M.n at [$>s] is identified with F^, S U', where 
P$ s is the orthogonal projector onto KerP^g (see lemma [3.3.2[) . The Hamiltonian differential 
equations for u are then equivalent to 

S-JL=-p* s U'. (1.22) 
or 

Given an initial value 'I's(O) = tpo E Sn, this equation has a unique solution r ^s(t) E Sn 
which satisfies the gauge condition (jl.7p . 

Main Theorem 1.7.2 (Adiabatic limit) Let be the smooth solution of (|1.20p . satisfying the 
gauge condition (|1.7|) . with smooth initial data ^(O), such that 

(i) lim^+oo 11^(0) — t/'oII'^ = 0; S or some smooth ipo E Sn, and 
(11) sup^ ||* M (0)|| Wa + ||^(0)|| Hl < K < 00. 
ITien i/iere exists r* > 0, independent of fi > 1, such that for s < 2, 

lim sup ||* M (r)-*s(r)|L =0 (1.23) 

where r 1— ► ^(t) £ 5^ is a curve m £/ie vortex space Sn, also satisfying (II. 7j) . which is the unique 
solution of (jl.22p wit/i initial data ^(O) = "00- r/ie projection onto the moduli space Mn : 

t h-> [*<j(t)] £7Wiv, 

is i/ie unique solution of the Hamiltonian system on (Sym N (T,),Q) associated to the Hamiltonian 
u defined in remark \1.7.1\ with initial value [4>q] £ Mn- 

This theorem in proved in employing a strategy which is explained in §1.81 following discussion 
of a very simple model problem. Some of the novel features which arise in the implementation of 
this strategy for (|1.18p are highlighted at the begininng of $2j 

Remark 1.7.3 (Related work) The approximation of the dynamical system (|1.18p by a dynam- 
ical system through a space of equilibria (in this case the self-dual vortices, which are the equilibria 
for A = 1) is referred to as an adiabatic limit or approximation. It was suggested in [30], following 
earlier conjectures of the same author on vortex and monopole dynamics in second order Lorentz 
invariant systems discussed in [31] . Proofs of the validity of the approximation in the case of second 
order dynamics were given in [45^146]: the strategy for the proof here, however, is different from that 
adopted in those references - see the discussion in §1.81 There has also been work on corresponding 
problems for cr-models, see |20|. 136] . A review of the analysis of adiabatic limit problems is given in 
[48] , mostly directed towards infinite dimensional natural Lagrangian systems of the type appearing 
in classical field theory. (Natural Lagrangian systems are those derivable from Lagrangians of the 
classical "kinetic energy minus potential energy" form). 
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Remark 1.7.4 (Implications for Chern-Simons vortex dynamics) Although it is not gen- 
erally possible to evaluate explicitly the Hamiltonian and symplectic form in the reduced system 
(jl.22p . it is possible to understand some basic features of the vortex dynamics in this model, see 
[271 l30t 139} l3l] . This work has been directed mostly to the case when the spatial domain is R 2 , so 
our theorem 11.7.21 does not imply the validity of the approximation (jl.22p in this case, see below. 
One general conclusion is that in the Chern-Simons model a force acting on the vortex produces 
motion at right angles to the direction of the force (in distinction to the behaviour in the relativistic 
case [31, 45J). Now it is known computationally (see |28} 131] and references therein), and in some 
special cases analytically (|47|). that the potential energy between two vortices depends on the 
distance between them, and is attractive for A < 1 and repulsive for A > 1. From this it can be 
deduced that two vortices will circle about one another, the direction of rotation depending upon 
whether A < 1 or A > 1. See [31] Section 7.13] for a discussion of these solutions in the R 2 case. 
Also in the same reference it is observed that (11.220 possesses another related type of solution: 
a rigidly rotating p-gon, with p vortices placed at the vertices of a regular p-gon. Many of the 
arguments and calculations leading to the conclusions about vortex dynamics can be carried out 
equally well with spatial domain the standard sphere £ = S 2 (|37j). even with explicit formulae in 
special limiting cases ([47]), in which case theorem II. 7. 21 imp lies rigorously the rotational behaviour 
for vortices described above. In future work results on the existence and stability of such periodic 
solutions for the full system (|1.5j) will be presented. 

It is to be hoped that some of these qualitative conclusions about vortex dynamics, (which are 
justified for (II. 5ft by the Main Theorem 1 1.7. 2 p would have a wider validity for Chern-Simons models 
of vortex dynamics, not necessarily close to any self-dual limit. There is some numerical evidence 
for this in related situations, for example the scattering of vortices in the relativistic abelian Higgs 
model is qualitatively similar for all values of the Higgs coupling constant, even though a rigorous 
analysis in which the vortices actually collide is only possible in the self-dual limit; see [45, 31J). On 
the other hand, the case of first order dynamics is in some ways numerically more problematic since 
it is not possible to produce any motion via choice of initial conditions (as can be done in the second 
order case), and it is necessary to have A deviate from the self-dual value 1, and quite substantially 
so in order to get motion which is easily computationally observable. A numerical study in [27] 
which compares the approximation (|1.22p with a computer simulation of (II. 5p finds that, in the 
case of spatial domain T, = R 2 , while the qualitative behaviour of two vortices is similar to that 
implied by (jl.22p for |A — 1| small, there are quantitative differences between the full dynamics and 
the adiabatic limit, which become quite marked as A moves away from the value 1. As the authors 
of [27] say, it is unclear to what extent some of these differences are genuine errors due to the neglect 
of radiation in the finite dimensional truncation (|1.22l) . as compared to being a numerical artefact; 
certainly some of the observed behaviour is consistent with energy being transferred into radiative 
modes, causing the vortices to spiral in towards one another in the attractive case ([271 Figure 6]). 
In any case, there is no issue with radiation when £ is a compact spatial domain, in which case 
theorem 11.7.21 does imply the validity of the approximation (11.22H for sufficiently small | A — 1 1 , and 
it seems reasonable to expect that in this case the dynamical bevaviour predicted by our analysis 
(relating (jl.5p to (jl.22p for small |A — 1|) is at least qualitatively relevant to the applications in the 
theoretical physics literature. 
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1.8 A simple model problem and discussion of methodology 



We consider here a simple two-dimensional example in order to exhibit as clearly as possible the 
phenomenon under study, and the strategy which will be employed in the proof of theorem 1 1.7.21 (It 
is the basic strategy taken in [40J for finite dimensional natural Lagrangian systems, here adapted to 
the case of infinite dimensions and to take advantage of the Bogomolny structure.) For real numbers 
(5 and (i^>l, we consider a linear first order Hamiltonian system for z(t) = (z (t), z 2 {r)) £ C 2 : 

Theorem 1.8.1 For each /j,^$>1, let r i— > Z^ir) € C 2 be the solution of 



i 1 = +p z 2 ) 

z 2 = i^z 1 + fiz 2 ), 



(1.24) 



with initial data satisfying |(Z*(0), Z 2 (0)) — (7, 0)| = 0(fi 1 ) as fi — > +oo, for some fixed 7 G C. 
Then 

lim wbx\ZJt) - (7e iT ,0)| = 0. (1.25) 

Remark 1.8.2 The system (|1.24p is Hamiltonian with the standard symplectic structure on C 2 
and with Hamiltonian function + U with V(z) = \z 2 z 2 and 

U(z) = -*V + + z 2 z l ). 

Thus V acts as a constraining potential for h — ► +00, forcing the solution onto the set S = C x {0} C 
C 2 where z 2 = 0. Projecting the system to S gives, formally, 



iz 1 + z 1 



0. 



(1.26) 



The theorem asserts that (jl.26p indeed governs the behaviour of the limit of appropriate sequences 
of solutions to (|1.24j) . 



Proof The solution with initial data z(0) = (z (0), z 2 (0)) is given by: 



z\r) 
At) 







(1-A_)e 



(l-A+)e 



iX - A z 1 {0) + (3(e iX+T -e iX ~ T 



z\0) 



1 



/5(A+ - A_) 



(l-A + )(l-A_)(e iA+T - t 



+ 



< x - T )z\o) 

iX-LT ( -\ \ \A\—T 



(1 - A+)e lA + T - (1 - A_)e* 



0(X+ - A_) 

Here the A± are the characteristic values of the system: 

4(^-/? 2 )V 



z 2 (0) 



1± 1 



A± = -(l + /u) 



which satisfy, by the binomial expansion, 

|A+ - n\ = 0(1), |A_ 



1 



as ^ — > 00. From this, and the fact that \± € M for large fi so that |e* A±r | = 1, the behaviour in 
(I1.25P follows for the solutions Z^{t) with initial data as described. □ 
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Remark 1.8.3 In this example the exact solutions indicate that while Z 2 — > 0, the time derivatives 
Z 2 are bounded, but cannot generally be expected to have limit zero. 

In the absence of explicit formulae for Z^{r), it is still possible to prove results like theorem 
IDTT1 either 

(i) by explicit perturbative construction of solutions to the full system, using solutions of the 
restricted system as a starting point, or 

(ii) by obtaining uniform bounds for the Z^{j) which allow the extraction of convergent subse- 
quences, and then identifying the unique limit of all such subsequences as the corresponding 
solution of the restricted system with Hamiltonian U\ s . 

In the present article we will adopt the second strategy in our proof of theorem 11,7.21 (although 
it would be possible to use the first strategy, as in [IS]). To make the structure of the proof 
transparent, it is useful to consider in some detail how to execute the second strategy to prove a 
variant of theorem 11.8.11 

Theorem 1.8.4 (Weaker version of theorem 11.8. ip In the situation of \1.8.1\ 

lim max |Z„(r) - ( 7 e ir ,0)| = 0, (1.27) 

/i— >+oo a<r<b 

for every bounded interval [a, b] C M. 

Remark 1.8.5 Although weaker than theorem ll.8.1l the proof of theorem ll.8.4l that we give gener- 
alizes to the infinite dimensional problem (jl.5p . (|1.18p . in which the explicit solutions corresponding 
to those used in the proof of theorem 11.8.11 are of course not available. 

Proof 

• Differentiation of the equations (jl.24jl in time gives the identical system £ = i. Use the energy 
identity: 

/xV(C(r)) + U(((t)) = AtV(C(0)) + *7(C(0)), 

together with the identical estimate for z(t), to deduce (using Cauchy-Schwarz) that the 
solutions Z^ of theorem 11.8.11 satisfy |Z^(r)| + |Z^(r)| < C, with C independent of /US>1. 

• By the previous item, deduce that the family of functions r i— > Z^{t) is uniformly (in /^3>1) 
bounded and equicontinuous, and so the Arzela-Ascoli theorem implies subsequential conver- 
gence Zfj,. — > Z in C(I) for any bounded interval / CM. 

• The energy estimate implies that, for large \i there exists C > 0, independent of fi, such that 
fj,Z 2 Z 2 < C. It follows that Z 2 — ► along any convergent subsequence. Now consider the 
integrated form of the first equation of (|l,24p (i.e. project the system onto S = C x {0} C C 2 
where z 2 = 0). Taking the limit — > oo, it follows that the limit Z = (Z 1 ,^ 2 ) of any 
convergent subsequence satisfies Z 1 ^) = i J Q r Z 1 (r')dr' and Z 1 (0) = 7. This integral equation 
has unique solution Z 1 {t) = je lT , and hence the C; oc limit of any convergent subsequence is 
(7e* r , 0). It follows that Z^ converges to this limit in Ci oc without restriction to subsequences. 
This proves theorem 11.8.41 (In view of remark 11.8.31 we should not expect this convergence 
to be in C} oc .) □ 
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The general situation to which theorem (jl.8.4p . and its proof, potentially generalize is the 
following: on a phase space 7i we consider the integral curves Z^t) for a Hamiltonian fiV + U for 
large jj, ("the full system"). Under the assumption that S = {z £ Ti : min V = V(z)} is a symplectic 
submanifold of TC, we can consider the "restricted system" on S determined by the Hamiltonian 
U\ s , and try to prove that this Hamiltonian system can be used to describe the limiting behaviour 
of Z^(t) as fi — > +00. An infinite dimensional example of this situation is provided by the Chern- 
Simons-Schrodinger system ()1 . 18|) : in the next section we will provide a proof of theorem 11.7.21 
employing the same strategy to that used in the proof of theorem 11.8.41 just given. 

2 Uniform bounds and proof of the main theorem 

In this section we prove our main result, theorem 11.7.21 along the lines suggested by the discussion 
of the simple model problem in the last section. The crucial stage is the proof of the main estimate, 
theorem l2.3-H which asserts the existence of a time interval, independent of fj,, on which the solution 
if) = (A, $) is uniformly bounded in Ti.2, and its time derivative is uniformly bounded in H 1 as 
fi — > +00. Given this bound, theorem 11.7.21 can be deduced using a variant of the Lions- Aubin 
lemma, and a careful analysis of the /U — ► +00 limit of (jl.18)) . Before obtaining the uniform bound, 
we collect some identities used in the proof. Some more specialized identities related to the self-dual 
structure are collected separately in $3j and referred to as needed. Specifically, we draw the reader's 
attention to the following two uses made of these more specialized identities: 

(i) Differentiation in time gives rise to an equation (|2.3ip for £ = tfj in which the dominant term 
(as fi — ► +00) involves the Hessian of V defined in (|2.40p . It is shown in $3] that this 
operator takes the special form 

L^=V* i p i , + 0{\B\), (2.28) 

with T>^ complex linear (see (|3.62p ). and B as in remark ["1.6. 3 1 Observing that the L 2 norm is 
exactly preserved for equations of the form J£ = V^V^^, it is easy to believe that the stated 
structure of is useful in the derivation of /i-independent bounds for (|2.3ip . (for initial data 
as in the theorem) ; this indeed turns out to be the case - see the proof of theorem 12,3. 11 

(ii) After obtaining a convergent subsequence of solutions of (|1.20p it is necessary to take the 
limit of the equation itself along the subsequence fj, = Hj — ► +00. For this purpose it is very 
convenient to be able to eradicate the term fiV' on the right hand side, since this is clearly 
hard to control for large \i: this can be done by applying a projection operator whose 
existence close to the set of self-dual vortices is assured by the Bogomolny structure: see 
lemmas f3 . 3 . 1 1 and 13.3.21 (In geometrical terms there is a foliation of the phase space TL2, and 
the range of is the tangent space to the leaves of this foliation, after dividing out by the 
action of the gauge group using (II. 7j) .) 

Although our final conclusions are in terms of the standard Sobolev norms based on the fixed 
connection V, it will be convenient to obtain bounds for the corresponding Sobolev norms defined 
at each fixed time with respect to the connection D = V — iA., see (|A.2p . These can be related to 
the standard norms by (!A.3P - (lA.5p . 
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2.1 The evolution equations and associated identities 

In addition to the rescaled equation (|1.20p for ip = (A, <&): 




= ixV + U' + J(dA ,iA $), 



we will use the differentiated equation for ( = tp = To write this down we need the linearization 
of the operator V (?/>), i.e. the second order linear differential operator Lw, obtained by differentiation 
of the map V' l— ► V(V0 : 

= DV'ty), 

or equivalently, (£,L^£) L 2 = ^V(V> + sC) U=o- Explicitly, with ( = (A, 6), we have 

(C, LyCiv = J (\dk\ 2 + \D$>\ 2 + |$| 2 | A| 2 - 2(D$,iA$) - 2(D6, iA$> (2.29) 

+ ($,$)2-i(l-|$|2)|$|2^ /ig . 



Remark 2.1.1 There is a slightly simpler version of this formula, given in (|2.40p below, when £ is 
restricted by the gauge condition (jl.7p . Furthermore in $3] it is shown that the self-dual structure 
provides a useful way of rewriting this formula as in (|2.28p . in terms of the complex structure 
defined in (jl.2ip . and using the complex one- form adz, where a = in place of the real 

one-form A\dx l + A^dx 1 , see (|3.60p . Since this is used only at one point in the proof - in lemma 
12.3.81 - this formulation is presented separately in £j3j and referred to only as needed. 

The linearization of U' is the linear operator Kj, = DU'(ip), given by 

K$ = (A, 6) ^ (o, ~(1 - |$| 2 )6 + , (2.30) 

with a defined in (|1.17p . Given these definitions, the chain rule implies that, if if) is a smooth 
solution of (|1.20p . then £(t) = V>(t) solves 

dC d 
$-±= f jiL i ,( + K i) t + 3-Q-(dAo,iAo$). (2.31) 

We also need identities for the evolution of the Bogomolny operator B defined in remark 11.6.31 
and discussed in more detail in The first component is preserved 

JL((£ _ 1(1 _ |$|2)) = e -2p( dl A 2 - d 2 A 1 ) + 6) = 0, (2.32) 

as a consequence of (|1.18p . We will require that the initial data are such that B — 1(1 — |$| 2 ) = 
initially, and hence for all times. The second component of the Bogomolny operator B will be 
denoted 

r/ = <9 A $ = ipi+iD^, (2.33) 
(see 53|), and we have the following identity: 

i(d T - iA )v = M-4<9 A (e- 2 ^ A 7?) + l$N " ^9 A ((1 - |$| 2 )$) . (2.34) 
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(To verify this identity: substitute A A $ = 4e~ 2p dAdA® - B$ into the third line of (|1.18j) and 
then apply 8a to the resulting equation and use the identity (E% + iE 2 )<& = — 2^|<£| 2 <9a < £ which 
follows from the first two lines of (|1.18p .) 
Of course, the energy 

£(t) = /uV(V(r)) + U(iP(t)) = £q>0 (2.35) 
is independent of time r for regular solutions, as is the L 2 norm 

||$(r)||L2=L>0. (2.36) 

2.2 Choice of gauge condition and related estimates 

The divergence of E can be calculated to be: 

div£ = e- 2p (d 1 E 1 + d 2 E 2 ) 

= J(-AB - e- 2p d 1 (i$,D 2 $) + e~ 2p d 2 (i$, D&) 



= ^e~ 2p \ V \ 2 ) + (^ - iA )<f>) - aB\<Z>\ 2 . 

In the last line we have used B = \(l - |$| 2 ), so that AB = -($, A A $) - e" 2p (|L»i$| 2 + \D 2 $\ 2 ), 
the equation for <3? and the definition of 77 in (|2.33p . Under the gauge condition (jl.7p we get the 
following equation for Aq: 

(-A + |$| 2 )A) = 4fie- 2p \r]\ 2 - aB\<S>\ 2 . (2.37) 



Lemma 2.2.1 (Estimates for Aq) Assume r i— > ip(r) = (A(r),$(r)), is a smooth solution, of 
(|1.20p which satisfies the gauge condition (|1.7p . f|2.35j) and (|2.36p . Then for all r < oo, there exists 
co(£o,L,r) > such that, 

\\A (r)\\ Lr < c (£ ,L,r) (2.38) 
and there exists cq (£q , L) > smc/i £/iai 

Po(r)||^< co(fo,i)(l + Ml|SA*(T)||r«). (2-39) 

Remark 2.2.2 This shows that in the original system (before rescaling) the time component of 
the potential Aq is 0(\X — 1|) in the gauge defined by (|1.7p . 

Proof The crucial point here is the \x independence of the bounds. The second inequality follows 
from standard elliptic theory once the first is established. By ()2.37j) it is possible to write Aq = 
Aq + Aq where (—A + l^l 2 )^^" = 4fj,e~ 2p \rj\ 2 , so that Aq > by the maximum principle, and 
(—A + |$| 2 )Ao = — <7-B|<3?| 2 . The bounds stated in the lemma will follow by the triangle inequality 
once they are proved for Aq , since they are immediate for Aq. Now integrating the equation for 
Aq implies that |||#| 2 Aq = J* s |<3?| 2 ^4||(i/i g < C(£q,L) since Aq > 0; this bound is independent 
of \i S> 1 on account of (12.35H . The standard elliptic theory for —Au = f G L 1 now gives the U 
estimates for At and hence the lemma. □ 
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Lemma 2.2.3 (Estimates for A) Let £ = (A,$) satisfy the gauge condition (jl.7p . as well as 
the linearized constraint equation (I2,32|) , Then there exists a constant c\ > such that ||A||#i < 
ci || $$11^2, and more generally, for any 1 < p < do, i/iere existe a constant ci(p) > suc/i i/iai 
||A|| w i, P < cx||$$||lp. in particular these estimates hold for a smooth solution, r i— > V( T ) = 
(A(r), $(t)), o/ (jl.20p which satisfies the gauge condition (|1.7|) . 

Proof These are the standard estimates for the Hodge system, proved by using the Hodge decom- 
position to reduce to the Calderon-Zygmund estimate for the Laplacian. □ 
On the subspace of ( = (A, satisfying the gauge condition (II, 7jl . the operator L^ has a 
simpler form: = L^(, where is the operator defined by 

(C,I^0l2 = J (W| 2 + |div A| 2 + |D$| 2 + |$| 2 (|A| 2 + |6| 2 ) (2.40) 



-4(D$,iAi)--(l-|^| 2 )|$| 2 )^ 9 . 



1 

2 ( 

Lemma 2.2.4 (The Hessian) Let ip = (A, <&) 6e smooth. Then the second order differential 
operator L^ is a self-adjoint operator with domain H 2 , and there exist numbers C2,c% such that 



<C,AK>^>c 2 ||Cll^i -c 3 ||Cll2 2 . 

The numbers C2,c^ depend only on the numbers L and So, defined as in (|2.35|) . (|2.36|) . 
Proof First of all, observe that 

j (\dk\ 2 + |div A| 2 + |D$| 2 + |$| 2 (|A| 2 + |#| 2 )W > c(S ,L)\\(A,^)f H 



This can be proved by a straightforward contradiction argument that is very similar to the proof 
of lemma 13.2.21 given below, so the details will be omitted. Next, to deduce the stated result, just 
bound the final two terms in (|2.40p using the Holder inequality with 1 = h + \ + \, the interpolation 
inequality in lemma lA. 91 and Cauchy-Schwarz. □ 

Corollary 2.2.5 Assume given a smooth solution, r V'( r ) = (A(t), $(t)), °f fjl .201) which 
satisfies the gauge condition (jl.7|) . (|2.35|) and (|2.36p . Then the quantity 

£i{r) = \&{t), {Lip + ^ 1 K^)C(t)) L 2, (2.41) 
where if) = iP(t), satisfies for \x > 1 

£x > c 4 \\c\\h A -<*K\\b 

with C4,c$ depending only on £q,L. 
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2.3 The main estimate 

We say that a smooth solution, r i— > V( r ) = C<4( r )> ^K t ))j 01 (|l-20p satisfies conditions (AE) and 
(AI), if the following conditions hold: 

(AE) There exists positive numbers £o,L such that ||$(r)||x,2 = L and £ (r) = £o> for all 
times r £ R, where £ (r) is the energy (|2.35 j) . (Recall that both these quantities 
are independent of r.) 

(AI) The initial data are such that || ^(0)||^ 2 + 11^(0) ll^i < K < oo. (Recall the definition 
of the norms in (ll.lOp ). 

Theorem 2.3.1 For jj, > 1 let r i— ► ?/>(t) 6e a smooth solution of (ll.20p satisfying conditions (AE) 
and (AI), for some fixed numbers K,L,£q. There exist numbers r* > anc? M* > ; independent 
of fj,, such that 

(<KT),JU/>(T))_ m < M*. (2.42) 



max 

|t|<t. 



H2XH 1 



Beginning of proof of theorem \2.3.1\ By time reversal invariance it is sufficient to prove the bound 
for < r < r*, for some r* > independent of //. Let 

C(r) = ^(r)=^(r). 
For any M > ||C(0)IIl 2 there exists a time T(M, n) > such that 

sup ||C(r)|Ua < M. (2.43) 
0<t<T{M,h) 

We will prove that there exist positive numbers M#, r*, independent of /u, such that T(M*, //) > t*, 
and hence sup 0<T<Tt ||C( r )llz, 2 < Af*. The proof proceeds by obtaining a series of //-independent 
bounds, predicated upon (I2.43|) . which imply boundedness of (ip(T),ip(r)) in the Hilbert space H2 
defined in (jl.lOp for < r < r*. These bounds are now stated in a sequence of lemmas, all of 
which refer to a smooth solution of (ll.20p , (II. 7p which verifies ( AE) , ( AI) and (|2.43p for all r under 
consideration. 

Lemma 2.3.2 (Estimate for $ in H 2 ) There exists C\ = C\(£q,L) > 0, independent of pi, such 
that 

||$(r)|| Hi <C 1 (l + ||C(r)|U 2 )<C 1 (l + M). 
Proof Using the third equation of (I1.18P for $, we bound 

||A A *||ra < PIU2 + \\A $\\ L 2 + - |$| 2 )|| L 2. 

Now, by lemma IA.2.21 we can bound HVaVa^Hl 2 < H^a^Il 2 + c ( i ?o)||Va < ^||l 4 ) an d hence, 
by lemma IA.9I and Cauchy-Schwarz: || VaVa^Hl 2 < 2||Aa^ > ||l2 + c(£q, L). Therefore, using also 
lemma [272TT1 we deduce the bound ||<&(i)|| H 2 < c(l+ ||C(i~)IIl 2 ) < c(l + M), for some c = c(£ , L) > 
0, and the result follows. □ 



Corollary 2.3.3 3C 2 = C 2 (£ ,L) > such that, ||*(t)||x«» < C 2 (l + ^/ln(l + M)). 
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Proof This follows from lemma IA. Ill and the previous lemma. □ 
Lemma 2.3.4 (Energy estimate for ( = ip) There is a constant Cs(Sq,L) > such that, 

d£i 



(It 



< c 3 (i + ii*iiioo)iicn 2 m + c 3 |icii b L2 + c 3 |icn 4 L2 . 



(2.44) 



where E\ is the quantity defined in (j2.41|) . 

Proof Compute 4z£i, substitute from (|2.3ip . and use the observation that 

(K,(dA ,i$A )) L 2 = o, 

by the constraint equation B = 1(1 — in (jl.5p . to obtain 



(2.45) 



-± = (i$,iA <l>) L 2 + -{Q,[—,L i> + l r 1 K^\Q L 2. 

To handle the second term, we make use of the following bounds (written schematically, i.e. sup- 
pressing indices and inner products which play no role): 

II<K 3 IIli < HSIU-IICIbllCII^ < cll^lU-HCIIiallCllffi 
||*ava*||li < ||v A 4>|| L2 ||A|| L 4||6|| L 4 ^cll^ll^Hcil^llCIlK 2 
||<i> 2 vA|| L1 < ||cili4||VA|| L 2 <c||$||^||C||i 2 ||CII H i- 

All of these bounds follow directly from Holder's inequality, the interpolation inequality in lemma 
IA.2.11 lemma 12.2.31 and the bound 



II 4 



+ ||A||#1 < c||<l?||_L°c ||$|| L 2. 



It then follows, by inspection of the formulae for L^,,K^, in (|2.29p and (|2.30p . that the second term 
in can be bounded by a sum of terms of this type, and hence: 

<c, [§- , l^ +/ ,-%]c> l2 <c(i + pii 2 LOO )iicii^+c|iciii 2 + c |icii| 2 . 



Also, we can bound 



where r > 1 and 1/r + 1/r' = 1. Combining these with lemma [2.2.14 we obtain (|2.44p . completing 
the proof of the lemma. □ 

Corollary 2.3.5 There is a constant C 4 = C A {£ ,K,L,M) > such that, \\({t)\\ h i < C 4 (l + t), 
for all times r £ [0, T{M , fj,)] . 

Lemma 2.3.6 (Estimate for ij = c^a^) There exists C§ = Cs{£q) > such that, at each time r, 

HMIirt < C(\m Hl + \\A\\h + Ml-). (2.46) 
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Proof From the equation (|2.34[) for rj, and using the interpolation inequality in lemma IA.91 the 
elliptic term 

C {A ^^{-Ad A {e- 2 Pd A r ] ) + \^\ 2 r 1 ) 
satisfies, for some c = c{8q) > 0, 

H|£(A,*)»7llLa < 11*11*} + ll$IM|A|| L 2 + cP || L 4(l + IMI^ 2 ) + c||$|||oc. (2.47) 

We next see that (12.46P follows from the usual elliptic regularity estimate. Firstly, observe that 
associated to the operator £(a,$) is the quadratic form 

Q(A,*)fa) = (V,£(A,*)V)lHS)= [ (4|5 A r ? | 2 e - 4 " + |$| 2 |r ? | 2 e -^) dfi g , 

which is bounded below by cllr/ll 2 ! where c = c(Sq,L) > by lemma [3.2.21 It follows that 

A 

ll^llifi < c II^(a,«i>) ? ?IIl 2 i a result which can be strengthened by the following 
Claim: ||VaVa??||l 2 — c ll^(A,$) r ?llL 2 where c = c(£q,L) > 0. 
By the Garding inequality 

IIVaVa^Hl 2 < \\£(A,>s>)V\\lz + c(.£o,L)(\\Vav\\l* + \\v\\h\)- 

Finally, using the interpolation inequality (1A.9P and the Cauchy-Schwarz inequality, we deduce the 
inequality claimed. 

□ 

Corollary 2.3.7 There is a constant C 6 = Cq(£o, K, L, M) > such that, h\\8a${t)\\l°° < C 6 (l + 
r). 

Lemma 2.3.8 (Closing the argument: estimate for ( in L 2 ) There is a constant Ct(£q,L,M) 
such that C = jpf satisfies 

HCOOIlia < ||C(0)|| 2 2 e C7 i'o T (II^A*(s)|| L oo+||$(s)||2^)d^ 
Proof Compute, using (|2.3ip . that 

^HC(r)||| 2 = 2(JC,(^+^)C) 

since (by the gauge condition) (£, (dAo,i$Ao)) L 2 = 0, and ((, (0, iAo&)) L 2 = (using = 

pointwise). By corollary 13 . 2 . 1 1 and the formula for K^, there exists Cj = Cj{£q,L) > such that 

<CV(Ml|^*(r)|| i - + ||*(r)||ia.)||C(T)||i a 

and so the stated inequality follows by the Gronwall lemma. □ 
Completion of proof of theorem \2.3.1\ The previous lemma allows us to validate the claim that 
(|2.43p . and thus all the bounds in lemmas I2.3.2tf27378l in fact hold on a [i- independent interval 
[0, t*], thus closing the argument. Indeed, by corollaries 12.3.31 and 12.3.71 we have ^^^^(t)!!^^ + 
||$(r)||ioo < C 8 (l + r) for some C 8 = C 8 (£ ,L,M). Now let r*,M* be such that 

\\C(0)\\he C7C * iT * +T * /2) < Ml 
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(This is always possible for M* > ||C(0)IIl 2 and r* small.) Then it follows that (|2.43j) holds with 
T(M%,fj,) > r*, and that the bounds given in lemma [2.3.21 through corollary 12.3.71 hold on the 
interval [0,r*]. To conclude, we explain how to derive the bounds in (|2.42p . For £ = tp we have 
boundedness of ||C( r )ll#^ by corollary 12.3.51 Integrating in r gives the bound for ||A||#i in (|2.42[> . 
Also the Kato and Sobolev inequalities (|_28j) give a bound for <j> in LP for 2 < p < oo. Together 
with the boundedness of ||$||z,°° this implies boundedness of ||A||jyi, P by lemma [2.2.31 Hence, 
integrating in r and applying Sobolev's inequality we deduce boundedness of [|A||£,«>. Putting 
all this information into (|A.3|) , (IA.4[) we can deduce, from lemma 12.3.21 and corollary 12.3.51 that 
(<&(t), <&(r)) is bounded in the (r-independent) norm H 2 x H 1 as claimed in (ET42D . □ 

2.4 Proof of theorem 11.7.21 

There are three stages to the proof: 

• Deduce, from the uniform bounds of theorem 12.3.11 and the compactness lemma 12.4.11 that 
for any sequence fj,j — ► +oo, there exists a subsequence along which the Vl/^. converge. 

• Identify the limit of these convergent subsequences. 

• Deduce, from the uniqueness of the limit just identified, that the Vl/^ do in fact converge as 
\i — ► +oo (without restriction to subsequences). 

The first stage of the proof depends upon the following version of the Lions- Aubin compactness 
lemma (see |29[ lemma 10.4]), which is proved by a modification of the standard proof of the usual 
Ascoli-Arzela theorem: 

Lemma 2.4.1 Assume that (V,h) is a smooth vector bundle with inner product, over a compact 
Riemannian manifold (£, g), which is endowed with a smooth unitary connection V and correspond- 
ing Sobolev norms \\ ■ \\h s on the space of sections defined as in 133)/ . Assume that l,s are positive 
numbers with I < s. Assume f n ( T ) is a sequence of smooth time- dependent sections of V which 
satisfy 

max(||/ n (T)||H3 + ||/n(T)M <C. 

\t\<t* 

Then there exists a subsequence {fn j }'jLi which converges to a limiting time- dependent section 
f € C([-r*,T*];iI s (y)) J in the sense that, max| T |< r>f \\(f n (r, ■ ) - f(r, ■ -> 0, for every r < s. 

Applying this we infer immediately the existence of a subsequence fij — > +oo along which the 
solutions typj,. = (A^',$^) converge to a limit ^s( T ) i n the sense that 

lim sup \\%At)-* s (t)\L =0, (2.48) 

for r < 2. It follows from corollary (|2.3.7p . that 

lim sup II^am^^Hl 00 = 0, 

and since the other Bogomolny equation B = ^(1 — l^] 2 ) is satisfied as a constraint, we deduce by 
theorem II. 6. 1\ that *&s( T ) £ ^N, i- e - the limit ^s( T ) 1S a self-dual vortex for each r G [— r*,r*]. In 
addition, by (|2.42p we have 

ll^( T i) - ^(t 2 )||hi < M*|n - t 2 \ 
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so that, by (|2.48j) . the limit *$>s W1 U also satisfy 

W^sin) - ^5(^2) ll^w < c\n -r 2 | 

for r' < 1, i.e. the limit is Lipschitz, and in particular lies in W 1,00 ([— r#, r*]; L 2 ). 

For the second stage, we need to identify the limiting curve r \—* *f>s( T ) £ Sn as that described 
in remark 11.7.11 It is clear, from the conditions on the initial data in the statement of theorem 
11.7.21 that ^5(0) = V'o £ $n, arid so it remains to deduce the ordinary differential equation ()1.22|) 
which then determines the curve completely. To do this it is necessary to take the limit of f|1.20j) : 

JJ^ =nV' + U' + J(dA%,iA^) (2.49) 

as jj, — > 00. The first term on the right hand side is the most evidently problematic. However, 
since the limiting motion is constrained to the vortex space Sn, it is only necessary to take a limit 
projected onto the tangent space T^ s Sn- To this end, it is actually most convenient to introduce 
P At (r) = P* M ( r ) the spectral projection operator onto Ker 2?^ ( r ) = KerP^ , jD $(j ( t j, discussed 
in lemma 13.3.21 By the final statement of lemma 13.3.21 and the convergence of \I/ ^. in (|2.48p . 
we know that P^(t) converge, in the 1? — ► L 2 operator norm, to the operator P^, s ( r ), which is 
the spectral projection operator onto KerP$ s ( r ) = KeiV* S/s ^'Dq, s ( T y (This latter operator is also 

the orthogonal I? projector onto the tangent space T^ s Sn (subject to the gauge condition (jl.7p ). 
Apply the operator P^t) to the equation (|1.2(jp . to obtain: 

since J(dAo,iAo<&u) and V'(\Pu) are both in the kernel of P^ by lemma [3.3,21 We can now identify 
the limit of the right hand side as S ( T )U' (^ s{ T )) at each r, and the convergence is strong in L 2 (S), 
by (|2.48p and the above mentioned convergence of P^ (r) . For the left hand side it is necessary to 
consider the limit of the derivatives Noting that these are bounded in e.g. L 2 {[— r#, r*]; L 2 (S)), 
we may assume (by restricting to a further subsequence if necessary) , the weak in I? subsequential 
convergence to a limit which is the weak time derivative of 

U'^^)i 2 ([-r.,r*];L2(S)) ~» \7 j ~~Q~~i L 2 ([-t,,t,]\L 2 (T,)) > 

for every / 6 L 2 ([— r*, t*]; L 2 (E)). Now to identify the limit along a convergent subsequence 
/ij — ► +00, consider the projection operator Fy, s r T \. Choosing /(r, •) = P$ s ( T )(/(r, •)), and using 
the symmetry of P^ s ( r ) this implies that 

+T * _ /" +r * 5*, 



cw f +T * aw, 

</,P* s (r)J-^ 1 >L»(E)rfr = y_ <P* S (T)/, J-gf->^(S)dT 

d^s /" +r * 

% s (r)f, J-^r)L 2 (E)^ T = / (/, 1P ^s(r)J-^r)L2(S)rf r ; 



for any / G L 2 ([— r*, r*]; L 2 (S)). On the other hand, by the above mentioned convergence of P^(t) 
to P^ S ( T ) and the bounded convergence theorem we have 



(P w (r)/,J^) L2(s) (ir-(P^ (T) /,J^) i2(E) 



dr -> 0, 
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on account of the bound (|2.42p . Therefore, we have in the limit: 




(2.51) 



for any / G L 2 {[— r*, r*]; L 2 (E)). But since the limit is known by the above to be in W 1,00 {[—t^,t^];L 2 ), 
it is differentiable (with respect to r, as a map into L 2 ) almost everywhere (the standard result 
extends to Hilbert space- valued functions, see, e.g., [H prop. 6.41]); the derivative lies in the tan- 
gent space T\£, 3 Sn, which is the range of the projector Fy, s r T \. Consequently (|2.51|) implies that 
r i — ^ ^s(t) is a solution of (jl.22p . with equality holding in L 2 for almost every r. But this in 
turn implies that r i— > ^(t) is actually continuously differentiable into L 2 , and we have a classical 
solution of (fL22l) . 

Finally for the third stage: we have now identified the limit as a solution of the limiting 
Hamiltonian system specified using remark Tl. 7. 11 Choosing smooth co-ordinates on Mn as in [37] 
we see that this is a smooth finite dimensional Hamiltonian system, and as such its solutions (for 
given initial data) are unique. Therefore all subsequences have the same limit, and so we can assert 
full convergence without resort to subsequences. 

3 Equations and identities related to the self-dual structure 

Notation change: In this section time does not appear at all, and so the boldface A for the 
spatial component is not used: i.e. in this section only, A refers to the spatial part of the 
connection, A = Aidx 1 + Aidx 2 . 

Ginzburg-Landau vortices are critical points of the static Ginzburg Landau energy functional 
Va = fj,v\(A,&)dfi g introduced following (11.131) . The coupling constant A > is central to the 
theory of critical points of the Ginzburg-Landau functional and the value A = 1 is special as in this 
case the functional admits the Bogomolny decomposition introduced in remark ["l.6.3[ This allows for 
a detailed understanding of the critical points not available for general values of A, and the theory 
of critical points for such general values is incomplete. (There is, however, a substantial literature 
on the asymptotic behaviour of critical points in the A —* +oo limit, starting with [6]; see [U] and 
references therein.) This decomposition of V = Vi has proved to be very useful not only for the 
analysis of critical points, but also for the associated time-dependent equations of vortex motion. 
For our purposes we need in particular to derive a special form for the operator associated to 
the Hessian of V, see (I3.6ip . 

3.1 Complex structure 

To discuss the Bogomolny structure in detail it is useful to use a complex formulation so we 
introduce the complex co-ordinate z = x l + ix 2 for the complex structure J on S. Using this, 
there is a decomposition of the complex 1-forms = O > © Q 0,1 into the ±i eigenspaces of J, see 
notation 11.3.11 Let QP(L) be the space of p-forms taking values in the bundle L: then for p = 1 
there is a similar decomposition, 

fi 1 (L) = O 1 '°(L)®O - 1 (Z). 
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Applying this decomposition to D$> £ Q l (L) we are led to introduce the operator D ' 1 given by 

D°> 1 $ = -((V 1 -iA 1 )+i(y 2 -iA 2 ))<f>dz = d A <f>dz. 
For real 1-forms A\dx l + A 2 dx 2 £ this decomposition reads 

Aidx 1 + A 2 dx 2 = adz + adz, 

where a = Al ~^ A ? ; a nd the map A t— > a (resp. A t—* a) is an R- linear isomorphism from to fi 1 ' 
(resp. f^ 0,1 ), and ||j4|| 2 2 = 4 J* aae" ~ 2p dfi g . With this a notation we can write 

<9a^ = ~m — ia<l>. 

3.2 The Hessian 

The Bogomolny decomposition amounts to the observation that, with A = 1, 

V{A^)=V 1 {A^)= l -fj K 4\B A ^\ 2 e- 2 P + (i? - | (1 - |3>| 2 )) 2 ) dfM g + irN 
where N = degL. If the following first order equations, called the Bogomolny equations, 

8a$ = 0, 

1 o (3-52) 

B--(l- |$| 2 ) = 
2 V 1 1 ; 

have solutions in a given class, they will automatically minimize V within that class. 
We introduce the nonlinear Bogomolny operator associated to this decomposition, 

b : n^en°(L) — >njk©fi 0,1 (£) 

(b - ±(1 - \$\ 2 ) , ,9 A 

Using the norm ?])||?2 = J*(|/3| 2 + ^e~ 2p \r]\ 2 )dpig induced from the metric on the target space, 
wc sec that V(A,$) = \\\B{A,<S>)\\ 2 L2 +ttN as in remark ll.6.3t see [8] . The derivative of B at 
ij) = {A, $) is the map DB$ : fijj fi°(L) — ► O - 1 ^) given by 

i-» (*«Li +<*,*), 5a* - (3.53) 

where a = ' 4l ~ 2 ^ 2 and d = j4l ' 2 t ^ 2 . Using this complex notation allows a simple unified formulation, 
which takes account of the gauge condition (jl.7p : this condition is the real part of 

4e- 2p da-i<£>® = 0, (3.54) 

while the imaginary part of this expression is just the condition *dA + ($, $) = 0, appearing in the 
linearized Bogomolny equations. This suggests the introduction of the operators 

^ : (n^ e o°(L)) — > (nj, e o^l)) 
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given by 

P^d, 6) = (4e" 2 ^d - A $ - *a$) ,„ 

(3.56) 



We use the real inner product associated to the L 2 norms induced from the metric as above, i.e.: 
\d, 6), {a', $')) L2 = J (Ae- 2p ma' + k4<5') dfi g on ft 1 - 
'{p, V), (/?', ??')) L2 = f ( W' + 4e- 2 ^ryr ? ') c*M s on fig, ^(L) . 



Integrating by parts we deduce that 



P^(d, d>), 09, r?)^ = ^(d, *), X$(/3, » 7 ) / ^ 
so that P^, is the I? adjoint of T>^ and 

VVDAa, $) = ( - 9(4e" 2 ' , 5d - z$<i>) - i$(^i + zd$) , (3.57) 



- \e~ 2p d A {d A ^ - ia*) - i$(4e~ 2p ad + i<M 
5(4e" 2p ad) + i(d A $)$ + |$| 2 d , -Ae- 2p d A d A <S> + |$| 2 i> + i4e~ 2p ad A ^ 
We compare this expression with the operator defined in (12.401) : 

: (n 1 ' ©^ ^)) — ► (n 1 '°©n°(L)) (3.58) 
which defines the Hessian of V on the subspace on which the gauge condition (jl.Tf) is satisfied, i.e., 

<V>, Z^) i2 = P 2 V^, V0 = ^l,=oV(^ + e^), (3.59) 

for ip = (A,&) satisfying (jl.7p . Using mixed real/complex notation for A/a, (|2.40p implies the 
following formula: 

= ( - 45(e" 2p 5d) + |<I>| 2 d - Pi$) + D 2 $) , 

- A A <i> - - (1 - 3|$| 2 )6 + 2ie~ 2p A ■ . (3.60) 



Calculate A ■ = 2ad A & + 2ad A $ and -(i<3>, + P2*) = i®d A ® - i$d A <&, from which 
it follows that 

(L-W= " ! ^ $ . _ _ j. (3.61) 

V * ^ ^ ^ (5- |(1 - |$| 2 )) $ + 4ie- 2 ' , daA < I> / 

(Incidentally, observing that 

B(i4 + i, $ + $) = B(A, $) + P^ + Ql^l 2 , 
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with ip = (A, <3?) satisfying (jl.7p . the identity (|3.6ip can also be read off from the quadratic part of 
the Taylor expansion for V(A + A, <I> + $): 

±(iP,L^) L2 = \\V^\\ 2 + (b{,I>) , Q|#| 2 ,-^\ 

= \\v*i>\h + J E (\( B - 5(1 - l*l 2 ))l<i>l 2 + ^(B A ^ -ia*>) 



using the inner product on O > © S1 U (L) defined above. 



Corollary 3.2.1 Let J denote the complex structure defined in (|1.21|) . There exists a number 
c > 0, independent ofip= (a, <3?) and ( = ip = (a, <t>) G fi 1,0 © smc/i that 



L 2 



Proo/ By (EOfll) |(JC,A/>0L2 ~ (A/JC^Cl^I < |BW0U«»|Cl£a- Now the complex structure $ 
written in complex notation, i.e. acting on Q 1,0 © £1®L), is given by J(d,<3?) = (— ia, Corre- 
spondingly, on fijjL © r2 0,1 (L) we introduce the complex structure 3'(/3,rf) = (i(3,—irj). Then, by 
observation 

2tyJC = -tfT>i,C (3-62) 

Therefore, writing w = V^Q, we have {V^Q,V^C) L i = (—fw,w) L 2 = by skew-symmetry, and 
the result follows. □ 

Lemma 3.2.2 Assume there are positive numbers L,£q such that |$|^2 = L, and V\{A, <3>) = £q 
and A > 0. Then the quadratic forms 

Q%((3) = j 4\d(3\ 2 e- 2p + \<S>\ 2 \f3\ 2 dfi g on QSlfcand 

Q {A ^) = J s 4e ~ ip \ d AV\ 2 + e- 2p m 2 \v\ 2 dfi g on ^(L) 

are strictly positive, and in fact bounded below by (respectively) C\\f3\\ 2 Hl and C\\r]\\ 2 Hl where C is a 
positive number depending only upon the numbers L,£q. 

Proof We will present the proof for the quadratic form Q(a,&)(v) as the other is similar but easier. 
Clearly Q( J 4,$)( r ?) > and in fact Q(a,$)(v) = if and only if n = on S (because if Baij = 
then n has isolated zeros (as in [28], sec. 3.5); if $77 = then n = since $ = a.e. contradicts 
J E |<I>| 2 = L > 0. Furthermore, we show that Q(A,$)( r ?) > C M 2 2 for a constant c; to be precise there 
exists c = c(L,£q) such that 

Q(A,&)(j]) — c > f° r an ?7 such that 1 1 77 1 1 ^2 = 1. (3.63) 

We will prove this by contradiction. First we obtain some bounds. By gauge invariance we are free to 
assume that the Coulomb gauge condition divj4 = holds. With this gauge condition, we have the 
bound 1 1 A \ \ H i < c{£q) and so A is bounded in every LP space. Now use ||<9f?||Lp < H^a^Ulp + II^IIlp 
to deduce that 

L <C(l + Q(A,*)fa)) 
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for every p < 2, by Holder's inequality. This in turn implies, by the LP estimate for the inhomo- 
geneous Cauchy-Riemann system, that r\ is bounded similarly in L 4 , and so since A is also we can 
bound drj in L 2 and hence r\ in H 1 . Finally, since A and r\ are bounded similarly in L 4 , this imples 
that ||f?||^i < Cil + Q^ty (rj)), with C depending only upon £q, L. To conclude, in Coulomb gauge 

the A, $,77 are all bounded in H 1 in terms of L,£q, Q(a,$>)(v)- 

The contradiction argument now starts: assume (|3.63p fails. Then, by the bounds just obtained 
and the Banach-Alaoglu and Rellich theorems, there is a sequence (A„,§ u ,r] u ) with 

\\A v \\ m + ||V$„|| £ 2 < K(£o,L), 

\\$u\\l 2 = L an d ||^||i,2 = 1) such that 

Q(A„,<3>„){Vv) ► 

A v — ► A weakly in H 1 

Q v — ► $ weakly in H 1 and strongly in LP for any p < 00 
r\ v — > 77 weakly in H 1 and strongly in LP . 

This implies that = L > 0, Q(a,$){v) = which implies as above that $ = a.e. and 

contradicts as above that |$|/,2 is constant. This leads to 

Q(A,0)(v) > c i\v\h where c\ = ci(L,£ ). 
Finally just apply the bound above for [|.D77||£a to improve this up to the H\ lower bound claimed. □ 

3.3 The Bogomolny foliation 

We introduce a foliation associated to the Bogomolny operator, which we regard as a map between 
the following Hilbert spaces: 

b ■. ^(n^en ^)) — > L 2 (n^®n°>\L)), 

(4*)- (b-±(1-\$\ 2 ) , 8 A <^ . 

With this choice of norms B is a smooth function. The next result shows that it is a submersion if 
the energy is close to the minimum value: 

Lemma 3.3.1 There exists 6* > such that \\5a&\\l 2 < @* implies that KerD^, = {0}, and 
KerDq, is 2N dimensional (where N = degL). 

Proof D*(j3,Tj) = is equivalent to 

- 9/3 - i$fj = 

- Ae~ 2p d A r] - i$/3 = 0. 

Apply the operations Ad to the first and Ad a to the second of these equations to deduce that 

- Ae- 2p dd(3 + |$| 2 /? - Aie- 2p d A <&fj = 

- Ad A {e~ 2p d A ri) + ~ i(d A $)/3 = 
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The first two terms of these two equations are respectively the Euler- Lagrange operators associated 
to the quadratic forms Q${(3) and Qa,${v) studied in the previous lemma. Then we get the 
estimates 

Q*(P) < c\d A $\ L 2 \P\ l a\ti\l* 
Qa,M <c\8 a *\v\P\l*\v\l* 

which implies the result, since Qg>((3) > c\(3\'j I1 and Qa,${v) ^ c \v\ 2 H i ■ 1=1 
The natural geometrical context for the results of this section will now be explained. Define 
0* = {(A,$) G H 1 ^^ © tt°(L)) : ||<9a$||l 2 < 0*} which is an open set containing {tp = (A,$) : 
B{i>) = 0} C fT^fijj, © n°(L)). Furthermore, the previous lemma implies that T>^ : (0, ■ © 
fl°(L)) — > (O^. © Q 0,1 (L)) is surjective for tp G O*. By the discussion in the paragraph preceding 
(|3.55p . this implies that DB^ : © Q,°(L) — * Oj™ © 0°' 1 (L) is also surjective for ip G O*, and 
hence the level sets of B form a foliation of 0* whose leaves have tangent space equal to Ker DB^, 
by P3 §3-5 and §4.4]. The intersection of this tangent space with SC^ = {(A,$) : (A,<&) satisfies 
(01) } is Ker Pvp. 

Lemma 3.3.2 Assume tp G (ft 1,0 © f^(L)) n °*- The operators V*,V^ defined in (|3T571) are 
self-adjoint operators on I? , with domain H 2 , with 2N -dimensional kernel equal to KerT>^, and: 

ii^;^ciu 2 + nciu 2 >ciiciiif 2 . (3.64) 

Let P^, 5e f/ie orthogonal spectral projector onto KerV^V^ = KerV^. Then F^(V'(ip)) = anc ^ 
P^(J(dx, *X*a»)) = f or an U smooth real valued function \- Finally, if also ipv) G (Q 1 ' ®fl^(L))r\ 
and supj ||w 2 < 00 an ^ limj^+oo ||^)Cj) — = 0, /or a// r < 2, i/ie corresponding 

projectors P^o) converge to P^, m L 2 — ► L 2 operator norm. 

Proof The first assertion and the bound (|3.64p follow from lemma 13.3.11 and standard elliptic 
theory. The next statement follows by noting that if n G KerX?^,, then differentiation of V(ip) = 
±f\B{ip)\ 2 dn g + 7rN yields 



<„,y'(^)) i2 = l 



V(i/j + sn) = (B(^),DB i ,(n)) L 2=0 



s=0 



since KerX?^, C Ker DB^ by the discussion preceding (|3.55p . Next, n G KevV^ implies that 
F^(3(dx,ix^fj.)) = since integration by parts reduces this to the fact that n solves the first 
component of DB^n = in (|3.53p . 

The final statement follows by [25l § IV.3], if it can be established that Tj = X>* W) I^y) converges 
to T = D^D^p in the generalized sense of Kato (see |25} §IV.2.6]), or equivalently in the norm 
resolvent sense: 

lim Wii + T)- 1 - (i + Tj)- l \\ L 2^ L 2 =0. (3.65) 

j->oo 

To verify this convergence, it is convenient first of all to verify it in Coulomb gauge. So let 
■ipU) = (A^\ <£>(■?)) = e lXi ■ ifjV) a nd 'tp = (A, $) = e lx ■ ip be gauge transforms (as defined following 
(ll.lip ). such that divAV' = = div A The assumed properties of ip^ ensure that sup ||xjll/i" 2 < 00 
and that lim \\xj — x\\H r = 0, Vr < 2 so that also ip^ — » ip in TC r for r < 2. Now observe that in 
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Coulomb gauge the formula (|3.57p does not involve any derivatives of the connection one-form A at 
all. From this it is then immediate by inspection that (writing Tj = D^^D^j), and T = T>^T>^,) 

||(f - I-KHz- < SMWh* < cSjiWClli? + II^CIU*) (3-66) 

where 5j — > as j — > +00. But this last fact implies (by |25} Theorems IV. 2. 24-25]) that Tj 
converges to T in the generalized sense, and hence in the resolvent sense: 

lim W + f)- 1 -(i + f j )- 1 \\ L ^ L 2 = 0. (3.67) 

This would establish the convergence of the corresponding spectral projectors in Coulomb gauge. 
To go back to the original tpj it is just necessary to make use of the following gauge invariance 
property: on £ = (&, the induced action of the gauge group is g • (d,$) = (a,g<&) for any S 1 
valued function g, and 

f (e ix • C) = e ix • (TC) , 

and similarly with Tj, Xj replaced by T, \. This gauge invariance property implies that (i + Tj)^ 1 = 
e~ lX] ' o(i+Tj)~ 1 oe' lXj and (i+T)^ 1 = e~ %x o{i+T)~ l oe LX , where by o we mean operator composition, 
and e lx is shorthand for the operator e* x » etc. Finally, using lim \\xj — x\\h t = 0, Vr < 2 we see 
that (I3.66P and (|3.67p imply (I3.65p . completing the proof. □ 



Appendix 

A.l Operators 

To describe in detail the Laplacian operators which appear in the text, we assume S to be covered 
by an atlas of charts U a on each of which is a local trivialisation of L determined by a choice of a 
local unitary frame. (A smooth section of L then corresponds to a family of smooth functions 
$ Q : U a -> C so that on U a n Up we have $ Q = e^Qp with e ie ^ :U a nUp ->■ S 1 smooth.) We 
assume given a smooth connection D = V — iA on L acting as a covariant derivative operator on 
sections of L. Working in such a chart, and suppressing the index a, the Laplacian on sections $ 
of L is given by: 

- A A $ = -±-Dj{jty/gD&) = -e- 2 "(AA*). (A.l) 

This satisfies (-Aa*,*')^ = ||I D ( $ + ^OI^Uo- 

Next we need the Laplacian on one-forms. Starting with A = Aidx 1 + ^dx 2 € £7™, the 

negative Laplacian is the Euler-Lagrange operator associated to the Dirichlet form ^ f(|divA| 2 + 
\dA\ 2 )dfi g (with the norms inside the integral determined by g in the standard way). Transferring 
to complex form a = \{A% — 1A2) £ this Dirichlet form is just 1(a) = 8 f e~ Ap dada d[i g . The 
corresponding negative Laplacian —A 1 ' is then defined by (— A 1,0 a, /3)jya = 4-I(a + ef3)\ e= o where 
we use the induced inner product O 1,0 as in This leads to the following formula for the negative 
Laplacian —A 1 ' on a G Q ' : 

-A^a = -4d(e- 2p da), 
which is precisely the operator appearing in ^3j Similarly, on ft 0,1 (L) the negative Laplacian is 

-A°/r ] = -4d A (e- 2 Pd A r l ), 

which is the operator in (|2.34p . 
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A. 2 Norms and inequalities 

We define the Sobolev norms denned with the covariant derivative D = Va = V — iA. (We write 
Va in place of D for emphasis here.) The first Sobolev norm is defined by 

= I (l*| 2 + |V A <i>| 2 ) d/v (A.2) 

A JE 

In the above integral the inner products are the standard ones induced from h and g. The higher 
norms H\, ... are defined similarly, as are the W2 P norms for integral k and any p £ [1, oo]. The 
IP norms of the higher covariant derivatives arising from the connections Va and V are related as 
expressed schematically in the following: 

||V$||lp < ||VA*||iP + c||A||ioo||$|| iP , (A.3) 

[|VV$[|lp < ||VaVa$||lp + c||A|| L oo||V a $||lp (A.4) 
+c(l + ||VA <f>\\ LP + ||A|||oo||$|| LP ), 

||VVV$||lp < ||VaV a Va$||lp + c||A|| LO o||V a Va^||lp (A.5) 
+c(l + ||VA|| L oo + ||A||!oo)||V a $||lp 

+c(l + ||V 2 A|| LlJ ||$|| L r + ||A||!oo||$||lp), 

where q~ l + r _1 = p^ 1 . 

We now collect together some inequalities from |13| . 
The system of equations 

B = f divA = 5 (A.6) 

(where as above div : Q 1 — ► fi° is minus the adjoint of d) is a first order elliptic system which can 
be solved for A subject to the condition on J fdfi g dictated by an integer N, the degree of L. It 
can be rewritten 

dA = (/ - b)dfi g divA = # (A.7) 

and solved via Hodge decomposition as long as the right hand sides have zero integral. There is a 
solution unique up to addition of harmonic 1-forms which satisfies ||A||^i, P < c p (l + ||/||lj> + H^Hif) 
for p < oo. 

Lemma A.2.1 (Covariant Sobolev and Gagliardo-Nirenberg inequalities) For 

(£,S) as above and for (A, $) G {H 1 x ff A )( S ) then V A $ € £ 4 (£) and 

||V a ^||l4 < c||V A $||^ (A.8) 
and also for all 1 < p < oo, H\ <^-> W^ p L°° continuously on S. Also 

||V a $||l4 < cIIVa^II^^IIVa^II^ 2 + IIVaVa^II^ 2 ) (a.9) 

where c depends only on (T,,g). 

Lemma A. 2. 2 (Covariant version of the Garding inequality) For ^ = (A, $) such that the 
norms on £ appearing below are finite we have 

||V a V a $||l2 < ||Aa$|| L 2 +c\\B\\ 1 /*\\V a $\\ L 2 +c||$||^||V a ^||^ 2 ||VB||^ 2 (A.10) 
where c is a number depending only on (£,<?). 
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Lemma A. 2. 3 (Covariant version of the Brezis-Gallouet inequality) If A G -^(E) and 
$ G £/ien 

ll*IU-(S) < c(l + \\$\\ H i A ^Hl + im&J) (A.ll) 
where c depends only on (EflO- 

A. 3 Global existence results and different choices of gauge 

In this section we will summarize the existence theory for (jl.5p from [5] and [13], and explain how 
theorem 11.4.11 can be deduced from it. Existence theory can be worked out using various gauge 
conditions, and a choice of gauge is usually made to facilitate the calculations. The simplest condi- 
tion for the statement of the theorem, which also is convenient if we wish to make the Hamiltonian 
structure manifest - see §1.51 is the temporal gauge condition Aq = 0; however, the regularity is 
stronger in Coulomb gauge divA = 0. We have the following statements. 

Theorem A. 3.1 (Global existence in temporal gauge) Given data <5(0) G if 2 (£) and A(0) G 
iT (E), there exists a global solution for the Cauchy problem for (II. 5p satisfying Aq = 0, with reg- 
ularity $ G C([0,oo);# 2 (£)) nC 1 ([0 ) oo);L 2 (E)) and A G C 1 ([0,oo); H 1 ^)) . Furthermore, it is 
the unique such solution satisfying Aq = and satisfies the estimate 

\m)\\iPV)<ce^\ 

for some positive constants c,a,(3 depending only on (Y,,g), the equations, and the initial data. 

This can be derived from theorem 1.1 in [13], by applying a gauge transformation to put the 
solution obtained there into temporal gauge. To be precise the cited result gives a global solution 
(cio, a, cf>) of the system (jl.5j) satisfying the parabolic gauge condition ao = diva, and the gauge 
invariant growth estimate 

U\\ H ^)(t) < ce ae0t . (A.12) 

The solution satisfies G C([0, do); # 2 (X)) n C 1 ([0, oo); L 2 (E)) , a G C([0, oo); H 1 ^)) and a G 
C([0, oo); L 2 (S)). Now define X e C l ([0, oo); L 2 (E)) by d tX + a Q = and x(0) = 0. Define ($, A) = 
(4>e ltx ,a + <ix) : this gives a solution to (II. 5h satisfying the properties asserted in theorem IA.3.11 
(Most of this can be read off immediately, except perhaps to verify that A G C 1 ([0, oo); i7 1 (E)) , 
but this follows from the first equation in (|1.5j) . using the fact that Ao = and the right hand side 
is continuous into L 2 .) 

An alternative approach to local existence is given in [5], where it is shown that, in Coulomb 
gauge, systems of the type (II. 5p can be put in the form of an abstract evolution equation to which 
Kato's theory ([26]) applies. This yields the existence of a local solution denoted (A',$') with $>' 
continuous into H 2 on a time interval of length determined by the H 2 norm of the initial data. But 
the estimate (|A.12p above is gauge invariant, and allows continuation of the local solution to provide 
a global solution in Coulomb gauge with regularity G C ([0, oo); H 2 (T,)) n C 1 ([0, oo); L 2 (E)) and 
A' G C([0,oo);# 3 (£)) n C 1 ([0, oo); H 1 (E)) satisfying the Coulomb gauge condition divA' = 0. 

Finally, we explain how to obtain theorem 11.4.11 from these results. Given a solution A',$' in 
Coulomb gauge, as just described, define x(tj x ) to be the solution of 

(-A + |$'| 2 )x = div A' - $') = $'), 
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with x(0, x) = 0. Then it is easy to verify that A = A' + dx,<3? = &e lx satisfies (|1.7p . Under 
the condition ||$(t)||^ 2 ^ s ^ = L > the solution exists and is unique at time t; this condition 
is natural because ||$(£)[|l2(£) is independent of time for solutions of (|1.5p . Now by the above 
mentioned Coulomb gauge regularity and the basic estimates for the Laplacian we deduce that 
X £ C 1 ([0, oo); H 2 ). This gives the global existence theorem in the gauge stated in theorem 11,4. 11 
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